Mendelian randomization analysis of 34,497 German Holstein cows to infer causal associations between milk production and health traits.

BACKGROUND: Claw diseases and mastitis represent the most important health issues in dairy cattle with a frequently mentioned connection to milk production. Although many studies have aimed at investigating this connection in more detail by estimating genetic correlations, they do not provide information about causality. An alternative is to carry out Mendelian randomization (MR) studies using genetic variants to investigate the effect of an exposure on an outcome trait mediated by genetic variants. No study has yet investigated the causal association of milk yield (MY) with health traits in dairy cattle. Hence, we performed a MR analysis of MY and seven health traits using imputed whole-genome sequence data from 34,497 German Holstein cows. We applied a method that uses summary statistics and removes horizontal pleiotropic variants (having an effect on both traits), which improves the power and unbiasedness of MR studies. In addition, genetic correlations between MY and each health trait were estimated to compare them with the estimates of causal effects that we expected. RESULTS: All genetic correlations between MY and each health trait were negative, ranging from - 0.303 (mastitis) to - 0.019 (digital dermatitis), which indicates a reduced health status as MY increases. The only non-significant correlation was between MY and digital dermatitis. In addition, each causal association was negative, ranging from - 0.131 (mastitis) to - 0.034 (laminitis), but the number of significant associations was reduced to five nominal and two experiment-wide significant results. The latter were between MY and mastitis and between MY and digital phlegmon. Horizontal pleiotropic variants were identified for mastitis, digital dermatitis and digital phlegmon. They were located within or nearby variants that were previously reported to have a horizontal pleiotropic effect, e.g., on milk production and somatic cell count. CONCLUSIONS: Our results confirm the known negative genetic connection between health traits and MY in dairy cattle. In addition, they provide new information about causality, which for example points to the negative energy balance mediating the connection between these traits. This knowledge helps to better understand whether the negative genetic correlation is based on pleiotropy, linkage between causal variants for both trait complexes, or indeed on a causal association.

Genomic dissection of the correlation between milk yield and various health traits using functional and evolutionary information about imputed sequence variants of 34,497 German Holstein cows.

BACKGROUND: Over the last decades, it was subject of many studies to investigate the genomic connection of milk production and health traits in dairy cattle. Thereby, incorporating functional information in genomic analyses has been shown to improve the understanding of biological and molecular mechanisms shaping complex traits and the accuracies of genomic prediction, especially in small populations and across-breed settings. Still, little is known about the contribution of different functional and evolutionary genome partitioning subsets to milk production and dairy health. Thus, we performed a uni- and a bivariate analysis of milk yield (MY) and eight health traits using a set of ~34,497 German Holstein cows with 50K chip genotypes and ~17 million imputed sequence variants divided into 27 subsets depending on their functional and evolutionary annotation. In the bivariate analysis, eight trait-combinations were observed that contrasted MY with each health trait. Two genomic relationship matrices (GRM) were included, one consisting of the 50K chip variants and one consisting of each set of subset variants, to obtain subset heritabilities and genetic correlations. In addition, 50K chip heritabilities and genetic correlations were estimated applying merely the 50K GRM. RESULTS: In general, 50K chip heritabilities were larger than the subset heritabilities. The largest heritabilities were found for MY, which was 0.4358 for the 50K and 0.2757 for the subset heritabilities. Whereas all 50K genetic correlations were negative, subset genetic correlations were both, positive and negative (ranging from -0.9324 between MY and mastitis to 0.6662 between MY and digital dermatitis). The subsets containing variants which were annotated as noncoding related, splice sites, untranslated regions, metabolic quantitative trait loci, and young variants ranked highest in terms of their contribution to the traits` genetic variance. We were able to show that linkage disequilibrium between subset variants and adjacent variants did not cause these subsets` high effect. CONCLUSION: Our results confirm the connection of milk production and health traits in dairy cattle via the animals` metabolic state. In addition, they highlight the potential of including functional information in genomic analyses, which helps to dissect the extent and direction of the observed traits` connection in more detail.

Genomic evidence for the suitability of Göttingen Minipigs with a rare seizure phenotype as a model for human epilepsy.

Epilepsy is a complex genetic disorder that affects about 2% of the global population. Although the frequency and severity of epileptic seizures can be reduced by a range of pharmacological interventions, there are no disease-modifying treatments for epilepsy. The development of new and more effective drugs is hindered by a lack of suitable animal models. Available rodent models may not recapitulate all key aspects of the disease. Spontaneous epileptic convulsions were observed in few Göttingen Minipigs (GMPs), which may provide a valuable alternative animal model for the characterisation of epilepsy-type diseases and for testing new treatments. We have characterised affected GMPs at the genome level and have taken advantage of primary fibroblast cultures to validate the functional impact of fixed genetic variants on the transcriptome level. We found numerous genes connected to calcium metabolism that have not been associated with epilepsy before, such as ADORA2B, CAMK1D, ITPKB, MCOLN2, MYLK, NFATC3, PDGFD, and PHKB. Our results have identified two transcription factor genes, EGR3 and HOXB6, as potential key regulators of CACNA1H, which was previously linked to epilepsy-type disorders in humans. Our findings provide the first set of conclusive results to support the use of affected subsets of GMPs as an alternative and more reliable model system to study human epilepsy. Further neurological and pharmacological validation of the suitability of GMPs as an epilepsy model is therefore warranted.

Transcriptome driven discovery of novel candidate genes for human neurological disorders in the telomer-to-telomer genome assembly era.

BACKGROUND: With the first complete draft of a human genome, the Telomere-to-Telomere Consortium unlocked previously concealed genomic regions for genetic analyses. These regions harbour nearly 2000 potential novel genes with unknown function. In order to uncover candidate genes associated with human neurological pathologies, a comparative transcriptome study using the T2T-CHM13 and the GRCh38 genome assemblies was conducted on previously published datasets for eight distinct human neurological disorders. RESULTS: The analysis of differential expression in RNA sequencing data led to the identification of 336 novel candidate genes linked to human neurological disorders. Additionally, it was revealed that, on average, 3.6% of the differentially expressed genes detected with the GRCh38 assembly may represent potential false positives. Among the noteworthy findings, two novel genes were discovered, one encoding a pore-structured protein and the other a highly ordered ß-strand-rich protein. These genes exhibited upregulation in multiple epilepsy datasets and hold promise as candidate genes potentially modulating the progression of the disease. Furthermore, an analysis of RNA derived from white matter lesions in multiple sclerosis patients indicated significant upregulation of 26 rRNA encoding genes. Additionally, putative pathology related genes were identified for Alzheimer's disease, amyotrophic lateral sclerosis, glioblastoma, glioma, and conditions resulting from the m.3242 A > G mtDNA mutation. CONCLUSION: The results presented here underline the potential of the T2T-CHM13 assembly in facilitating the discovery of candidate genes from transcriptome data in the context of human disorders. Moreover, the results demonstrate the value of remapping sequencing data to a superior genome assembly. Numerous potential pathology related genes, either as causative factors or related elements, have been unveiled, warranting further experimental validation.

Blood transcriptome analysis in a buck-ewe hybrid points towards an nuclear factor-kappa B lymphoproliferative autoimmune disorder.

Mammal hybridization is a speciation mechanism and an evolutionary driver. Goat-sheep, especially buck-ewe hybrids, are very rare with only one case reported in 2016, which is the subject of the work presented here. Blood transcriptome analysis revealed that the hybrid largely deviated from imprinting schemes previously described in sheep and other mammals. Furthermore, transcriptome regulation seems to differ from the parent transcriptomes, which is most likely a product of partially incompatible imprinting mechanisms from two closely related species. To gain a deeper understanding of hybridization in mammals we re-analyzed the RNA sequencing data of the buck-ewe hybrid and its parents. We found parent-of-origin-specific expression of genes that functionally clustered, which we explain with the Dobzhansky-Muller incompatibility (DMI) model. According to the DMI model, proteins which interact have a high probability of being barrier loci and hence are prone to monoallelic expression. We discovered enrichment of genes uniquely expressed by the buck-ewe hybrid, which implicate that it suffered from an NF-κB lymphoproliferative autoimmune disorder. Similar findings were reported in the F1 generation of hybrid mice. We propose that hybridization of two related species may lead to an autoimmune phenotype, due to immunoglobulin incompatibilities and incomplete silencing of barrier loci.

Disturbed circadian rhythm of locomotor activity of pullets is related to feather pecking in laying hens.

Various aspects of activity, such as spontaneous activity, explorative activity, activity in open-field tests, and hyperactivity syndrome have been explored as causal factors of feather pecking in laying hens, with no clear results. In all previous studies, mean values of activity over different time intervals were used as criteria. Incidental observation of alternated oviposition time in lines selected for high (HFP) and low feather pecking (LFP), supported by a recent study which showed differentially expressed genes related to the circadian clock in the same lines, led to the hypothesis that feather pecking may be related to a disturbed diurnal activity rhythm. Hence activity recordings of a previous generation of these lines have been reanalyzed. Data sets of a total of 682 pullets of 3 subsequent hatches of HFP, LFP, and an unselected control line (CONTR) were used. Locomotor activity was recorded in pullets housed in groups of mixed lines in a deep litter pen on 7 consecutive 13-h light phases, using a radio-frequency identification antenna system. The number of approaches to the antenna system was recorded as a measure of locomotor activity and analyzed using a generalized linear mixed model including hatch, line, time of day and the interactions of hatch × time of day and line × time of day as fixed effects. Significant effects were found for time and the interaction line × time of day but not for line. All lines showed a bimodal pattern of diurnal activity. The peak activity of the HFP in the morning was lower than that of the LFP and CONTR. In the afternoon peak all lines differed with the highest mean in the LFP followed by CONTR and HFP. The present results provide support for the hypothesis that a disturbed circadian clock plays a role in the development of feather pecking.

Structural variation and eQTL analysis in two experimental populations of chickens divergently selected for feather-pecking behavior.

Feather pecking (FP) is a damaging nonaggressive behavior in laying hens with a heritable component. Its occurrence has been linked to the immune system, the circadian clock, and foraging behavior. Furthermore, dysregulation of miRNA biogenesis, disturbance of the gamma-aminobutyric acid (GABAergic) system, as well as neurodevelopmental deficiencies are currently under debate as factors influencing the propensity for FP behavior. Past studies, which focused on the dissection of the genetic factors involved in FP, relied on single nucleotide polymorphisms (SNPs) and short insertions and deletions < 50 bp (InDels). These variant classes only represent a certain fraction of the genetic variation of an organism. Hence, we reanalyzed whole-genome sequencing data from two experimental populations, which have been divergently selected for FP behavior for over more than 15 generations, performed variant calling for structural variants (SVs) as well as tandem repeats (TRs), and jointly analyzed the data with SNPs and InDels. Genotype imputation and subsequent genome-wide association studies, in combination with expression quantitative trait loci analysis, led to the discovery of multiple variants influencing the GABAergic system. These include a significantly associated TR downstream of the GABA receptor subunit beta-3 (GABRB3) gene, two microRNAs targeting several GABA receptor genes, and dystrophin (DMD), a direct regulator of GABA receptor clustering. Furthermore, we found the transcription factor ETV1 to be associated with the differential expression of 23 genes, which points toward a role of ETV1, together with SMAD4 and KLF14, in the disturbed neurodevelopment of high-feather pecking chickens.

Resistance to chicken amyloid arthropathy is associated with a dysfunctional mutation in serum amyloid A.

Chicken amyloid arthropathy is a debilitating disease with a major impact on animal welfare. Since the disease is triggered by bacterial infection, preventative treatment also contributes to the widespread overuse of antibiotics. Bacterial infection initiates an acute phase response including increased serum amyloid A (SAA) production by the liver. SAA accumulates at sites of infection and in particular in large joints of affected birds. Interestingly, white egg-laying chickens (WL) are resistant to the disease whilst brown egg-laying chickens (BL) are most affected. Disease susceptibility has an immunological basis but the possible contribution of underlying genetic risk factors is not understood. Using a whole genome sequencing approach, we discovered a novel variant in the SAA gene in WL, which is predicted to result in an arginine to serine substitution at position 90 (SAA.R90S). Surprisingly, when overexpressed in chicken hepatocellular carcinoma cells, SAA.R90S was expressed at a higher rate and secreted to a greater degree than the wild-type SAA protein. Moreover, RNASeq analysis showed that the R90S mutant exerted a differential effect on the expression of core transcription factors linked to cell fate determination and cell differentiation. Comparative analysis of gene expression in murine CD4 T-cells stimulated with IL-6/SAA, suggests that SAA.R90S might block an induced cell fate change toward pro-inflammatory T helper 17 cells, which are required for immunological protection against pathogenic bacteria during an acute phase response. Our results provide first mechanistic insights into the genetic resistance of WL to amyloid arthropathy and could be applied to commercial layer breeding programs to improve animal welfare and reduce the negative effects of the overuse of antibiotics.

eQTL analysis of laying hens divergently selected for feather pecking identifies KLF14 as a potential key regulator for this behavioral disorder.

Feather pecking in chickens is a damaging behavior, seriously impacting animal welfare and leading to economic losses. Feather pecking is a complex trait, which is partly under genetic control. Different hypotheses have been proposed to explain the etiology of feather pecking and notably, several studies have identified similarities between feather pecking and human mental disorders such as obsessive-compulsive disorder and schizophrenia. This study uses transcriptomic and phenotypic data from 167 chickens to map expression quantitative trait loci and to identify regulatory genes with a significant effect on this behavioral disorder using an association weight matrix approach. From 70 of the analyzed differentially expressed genes, 11,790 genome wide significantly associated variants were detected, of which 23 showed multiple associations (≥15). These were located in proximity to a number of genes, which are transcription regulators involved in chromatin binding, nucleic acid metabolism, protein translation and putative regulatory RNAs. The association weight matrix identified 36 genes and the two transcription factors: SP6 (synonym: KLF14) and ENSGALG00000042129 (synonym: CHTOP) as the most significant, with an enrichment of KLF14 binding sites being detectable in 40 differentially expressed genes. This indicates that differential expression between animals showing high and low levels of feather pecking was significantly associated with a genetic variant in proximity to KLF14. This multiallelic variant was located 652 bp downstream of KLF14 and is a deletion of 1-3 bp. We propose that a deletion downstream of the transcription factor KLF14 has a negative impact on the level of T cells in the developing brain of high feather pecking chickens, which leads to developmental and behavioral abnormalities. The lack of CD4 T cells and gamma-Aminobutyric acid (GABA) receptors are important factors for the increased propensity of laying hens to perform feather pecking. As such, KLF14 is a clear candidate regulator for the expression of genes involved in the pathogenic development. By further elucidating the regulatory pathways involved in feather pecking we hope to take significant steps forward in explaining and understanding other mental disorders, not just in chickens.

Structural variants and tandem repeats in the founder individuals of four F2 pig crosses and implications to F2 GWAS results.

BACKGROUND: Structural variants and tandem repeats are relevant sources of genomic variation that are not routinely analyzed in genome wide association studies mainly due to challenging identification and genotyping. Here, we profiled these variants via state-of-the-art strategies in the founder animals of four F2 pig crosses using whole-genome sequence data (20x coverage). The variants were compared at a founder level with the commonly screened SNPs and small indels. At the F2 level, we carried out an association study using imputed structural variants and tandem repeats with four growth and carcass traits followed by a comparison with a previously conducted SNPs and small indels based association study. RESULTS: A total of 13,201 high confidence structural variants and 103,730 polymorphic tandem repeats (with a repeat length of 2-20 bp) were profiled in the founders. We observed a moderate to high (r from 0.48 to 0.57) level of co-localization between SNPs or small indels and structural variants or tandem repeats. In the association step 56.56% of the significant variants were not in high LD with significantly associated SNPs and small indels identified for the same traits in the earlier study and thus presumably not tagged in case of a standard association study. For the four growth and carcass traits investigated, many of the already proposed candidate genes in our previous studies were confirmed and additional ones were identified. Interestingly, a common pattern on how structural variants or tandem repeats regulate the phenotypic traits emerged. Many of the significant variants were embedded or nearby long non-coding RNAs drawing attention to their functional importance. Through which specific mechanisms the identified long non-coding RNAs and their associated structural variants or tandem repeats contribute to quantitative trait variation will need further investigation. CONCLUSIONS: The current study provides insights into the characteristics of structural variants and tandem repeats and their role in association studies. A systematic incorporation of these variants into genome wide association studies is advised. While not of immediate interest for genomic prediction purposes, this will be particularly beneficial for elucidating biological mechanisms driving the complex trait variation.

Development and validation of a horse reference panel for genotype imputation.

BACKGROUND: Genotype imputation is a cost-effective method to generate sequence-level genotypes for a large number of animals. Its application can improve the power of genomic studies, provided that the accuracy of imputation is sufficiently high. The purpose of this study was to develop an optimal strategy for genotype imputation from genotyping array data to sequence level in German warmblood horses, and to investigate the effect of different factors on the accuracy of imputation. Publicly available whole-genome sequence data from 317 horses of 46 breeds was used to conduct the analyses. RESULTS: Depending on the size and composition of the reference panel, the accuracy of imputation from medium marker density (60K) to sequence level using the software Beagle 5.1 ranged from 0.64 to 0.70 for horse chromosome 3. Generally, imputation accuracy increased as the size of the reference panel increased, but if genetically distant individuals were included in the panel, the accuracy dropped. Imputation was most precise when using a reference panel of multiple but related breeds and the software Beagle 5.1, which outperformed the other two tested computer programs, Impute 5 and Minimac 4. Genome-wide imputation for this scenario resulted in a mean accuracy of 0.66. Stepwise imputation from 60K to 670K markers and subsequently to sequence level did not improve the accuracy of imputation. However, imputation from higher density (670K) was considerably more accurate (about 0.90) than from medium density. Likewise, imputation in genomic regions with a low marker coverage resulted in a reduced accuracy of imputation. CONCLUSIONS: The accuracy of imputation in horses was influenced by the size and composition of the reference panel, the marker density of the genotyping array, and the imputation software. Genotype imputation can be used to extend the limited amount of available sequence-level data from horses in order to boost the power of downstream analyses, such as genome-wide association studies, or the detection of embryonic lethal variants.

Genetic and genomic characterization followed by single-step genomic evaluation of withers height in German Warmblood horses.

Reliability of genomic predictions is influenced by the size and genetic composition of the reference population. For German Warmblood horses, compilation of a reference population has been enabled through the cooperation of five German breeding associations. In this study, preliminary data from this joint reference population were used to genetically and genomically characterize withers height and to apply single-step methodology for estimating genomic breeding values for withers height. Using data on 2113 mares and their genomic information considering about 62,000 single nucleotide polymorphisms (SNPs), analysis of the genomic relationship revealed substructures reflecting breed origin and different breeding goals of the contributing breeding associations. A genome-wide association study confirmed a known quantitative trait locus (QTL) for withers height on equine chromosome (ECA) 3 close to LCORL and identified a further significant peak on ECA 1. Using a single-step approach with a combined relationship matrix, the estimated heritability for withers height was 0.31 (SE = 0.08) and the corresponding genomic breeding values ranged from - 2.94 to 2.96 cm. A mean reliability of 0.38 was realized for these breeding values. The analyses of withers height showed that compiling a reference population across breeds is a suitable strategy for German Warmblood horses. The single-step method is an appealing approach for practical genomic prediction in horses, because not many genotypes are available yet and animals without genotypes can by this way directly contribute to the estimation system.

The light response in chickens divergently selected for feather pecking behavior reveals mechanistic insights towards psychiatric disorders.

BACKGROUND: Feather pecking is a serious behavioral disorder in chickens that has a considerable impact on animal welfare and poses an economic burden for poultry farming. To study the underlying genetics of feather pecking animals were divergently selected for feather pecking over 15 generations based on estimated breeding values for the behavior. METHODS AND RESULTS: By characterizing the transcriptomes of whole brains isolated from high and low feather pecking chickens in response to light stimulation we discovered a putative dysregulation of micro RNA processing caused by a lack of Dicer1. This results in a prominent downregulation of the GABRB2 gene and other GABA receptor transcripts, which might cause a constant high level of excitation in the brains of high feather pecking chickens. Moreover, our results point towards an increase in immune system-related transcripts that may be caused by higher interferon concentrations due to Dicer1 downregulation. CONCLUSION: Based on our results, we conclude that feather pecking in chickens and schizophrenia in humans have numerous common features. For instance, a Dicer1 dependent disruption of miRNA biogenesis and the lack of GABRB2 expression have been linked to schizophrenia pathogenesis. Furthermore, disturbed circadian rhythms and dysregulation of genes involved in the immune system are common features of both conditions.

GWAS Hits for Bilateral Convergent Strabismus with Exophthalmos in Holstein Cattle Using Imputed Sequence Level Genotypes.

Bilateral convergent strabismus with exophthalmos (BCSE) is a malformation of the eyes and is recognized as a mild but progressive disorder that affects cattle in the first two years of life. This most likely inherited disorder is rarely described in cattle resembling autosomal dominantly inherited forms of human progressive external ophthalmoplegia (PEO). In German Braunvieh cattle, two linked genome regions were found that could be responsible for the development and/or progression of BCSE. The goal of this study was to phenotypically characterize BCSE in Holstein cattle from Germany and Switzerland as well as to identify associated genome regions by GWAS. The clinicopathological phenotype of 52 BCSE-affected Holstein cattle was in accordance with the phenotype described in German Braunvieh cattle, but in addition, signs of degeneration and cellular infiltration in the eye muscles were found. By using imputed sequence level genotype data, three genome-wide significant GWAS hits were revealed on different chromosomes that were not detected by initial GWAS based on high density SNP array data highlighting the usefulness of this approach for mapping studies. The associated genome regions include the ABCC4 gene as well as markers adjacent to the NCOR2 and DNAJC3 genes all illustrating possible functional candidate genes. Our results challenge a monogenic mode of inheritance and indicate a more complex inheritance of BCSE in Holstein cattle. Furthermore, in comparison to previous results from German Braunvieh cattle, it illustrates an obvious genetic heterogeneity causing BSCE in cattle. Subsequent whole genome sequencing (WGS)-based analyses might elucidate pathogenic variants in the future.

Differentially Expressed Gene Patterns in Ascarid-Infected Chickens of Higher- or Lower-Performing Genotypes.

Here, we describe the first transcriptomic investigation of the peripheral blood of chickens exposed to Ascaridia galli and Heterakis gallinarum infections. We investigated differentially expressed gene (DEG) patterns in two chicken genotypes with either a higher (Lohmann Brown Plus, LB) or lower (Lohmann Dual, LD) laying performance level. The hens were experimentally coinfected with A. galli and H. gallinarum, and their worm burdens and infection parameters were determined six weeks post infection. Based on most representative infection parameters, the hens were clustered into lower- and higher-infection intensity classes. We identified a total of 78 DEGs contributing to infection-related phenotypic variation in the two genotypes. Our data showed significant upregulation of Guanylate Binding Protein 7 (GBP7) in LD hens, making it a promising candidate for tolerance to ascarid infections in chickens. Gene ontology analysis revealed higher transcriptome activity related to biological processes such as "response to external stimulus" in LB hens, implying a higher stress response in this genotype. In contrast, LD hens showed higher transcriptomic expression of genes related to ontology classes that are possibly associated with a higher tolerance to infections. These findings may help explain why lower-performing genotypes (i.e., LD) are less sensitive to infections in terms of maintaining their performance.

Multi-species transcriptome meta-analysis of the response to retinoic acid in vertebrates and comparative analysis of the effects of retinol and retinoic acid on gene expression in LMH cells.

BACKGROUND: Retinol (RO) and its active metabolite retinoic acid (RA) are major regulators of gene expression in vertebrates and influence various processes like organ development, cell differentiation, and immune response. To characterize a general transcriptomic response to RA-exposure in vertebrates, independent of species- and tissue-specific effects, four publicly available RNA-Seq datasets from Homo sapiens, Mus musculus, and Xenopus laevis were analyzed. To increase species and cell-type diversity we generated RNA-seq data with chicken hepatocellular carcinoma (LMH) cells. Additionally, we compared the response of LMH cells to RA and RO at different time points. RESULTS: By conducting a transcriptome meta-analysis, we identified three retinoic acid response core clusters (RARCCs) consisting of 27 interacting proteins, seven of which have not been associated with retinoids yet. Comparison of the transcriptional response of LMH cells to RO and RA exposure at different time points led to the identification of non-coding RNAs (ncRNAs) that are only differentially expressed (DE) during the early response. CONCLUSIONS: We propose that these RARCCs stand on top of a common regulatory RA hierarchy among vertebrates. Based on the protein sets included in these clusters we were able to identify an RA-response cluster, a control center type cluster, and a cluster that directs cell proliferation. Concerning the comparison of the cellular response to RA and RO we conclude that ncRNAs play an underestimated role in retinoid-mediated gene regulation.

Meta-analyses of genome wide association studies in lines of laying hens divergently selected for feather pecking using imputed sequence level genotypes.

BACKGROUND: Feather pecking (FP) is damaging behavior in laying hens leading to global economic losses in the layer industry and massive impairments of animal welfare. The objective of the study was to discover genetic variants and affected genes that lead to FP behavior. To achieve that we imputed low-density genotypes from two different populations of layers divergently selected for FP to sequence level by performing whole genome sequencing on founder and half-sib individuals. In order to decipher the genetic structure of FP, genome wide association studies and meta-analyses of two resource populations were carried out by focusing on the traits 'feather pecks delivered' (FPD) and the 'posterior probability of a hen to belong to the extreme feather pecking subgroup' (pEFP). RESULTS: In this meta-analysis, we discovered numerous genes that are affected by polymorphisms significantly associated with the trait FPD. Among them SPATS2L, ZEB2, KCHN8, and MRPL13 which have been previously connected to psychiatric disorders with the latter two being responsive to nicotine treatment. Gene set enrichment analysis revealed that phosphatidylinositol signaling is affected by genes identified in the GWAS and that the Golgi apparatus as well as brain structure may be involved in the development of a FP phenotype. Further, we were able to validate a previously discovered QTL for the trait pEFP on GGA1, which contains variants affecting NIPA1, KIAA1211L, AFF3, and TSGA10. CONCLUSIONS: We provide evidence for the involvement of numerous genes in the propensity to exhibit FP behavior that could aid in the selection against this unwanted trait. Furthermore, we identified variants that are involved in phosphatidylinositol signaling, Golgi metabolism and cell structure and therefore propose changes in brain structure to be an influential factor in FP, as already described in human neuropsychiatric disorders.

Analysis of the brain transcriptome in lines of laying hens divergently selected for feather pecking.

BACKGROUND: Feather pecking (FP) in laying hens reduces animal welfare and leads to economic losses for the layer industry. FP is considered a heritable condition that is influenced by dysregulation of neurotransmitter homeostasis, the gut microbiome, and the immune system. To identify genes and biological pathways responsible for FP behavior we compared the brain transcriptomes of 48 hens divergently selected for FP. In addition, we tested if high feather peckers (HFP) and low feather peckers (LFP) respond differently to light since light has been shown to trigger FP behavior. RESULTS: Of approximately 48 million reads/sample an average of 98.4% were mapped to the chicken genome (GRCg6a). We found 13,070 expressed genes in the analyzed brains of which 423 showed differential expression between HFP and LFP. Genes of uncertain function and non-coding RNAs were overrepresented among those transcripts. Functional analyses revealed the involvement of cholinergic signaling, postsynaptic activity, membrane channels, and the immune system. After the light stimulus, 28 genes were found to be differentially expressed. These included an interaction cluster of core components of the circadian clock. However, differences in the response to light between HFP and LFP were not detectable. CONCLUSIONS: Genes involved in cholinergic signaling, channel activity, synaptic transmission, and immune response were found to be involved in FP behavior. We propose a model in which the gut microbiota modulates the immune system, which in turn affects cholinergic signaling. This might have an influence on monoamine signaling with possible involvement of GABA or glutamate signaling.

Blood transcriptome analysis in a buck-ewe hybrid and its parents.

Examples of living sheep-goat hybrids are rare, mainly due to incorrect chromosome pairing, which is thought to be the main cause for species incompatibility. This case represents the first report of a buck-ewe hybrid and the first mammalian hybrid to be analyzed with next generation sequencing. The buck-ewe hybrid had an intermediate karyotype to the parental species, with 57 chromosomes. Analysis of the blood transcriptomes of the hybrid and both parents revealed that gene expression levels differed between the hybrid and its parents. This could be explained in part by age-dependent differences in gene expression. Contribution to the geep transcriptome was larger from the paternal, compared to the maternal, genome. Furthermore, imprinting patterns deviated considerably from what is known from other mammals. Potentially deleterious variants appeared to be compensated for by monoallelic expression of transcripts. Hence, the data imply that the buck-ewe hybrid compensated for the phylogenetic distance between the parental species by several mechanisms: adjustment of gene expression levels, adaptation to imprinting incompatibilities, and selective monoallelic expression of advantageous transcripts. This study offers a unique opportunity to gain insights into the transcriptome biology and regulation of a hybrid mammal.

GWAS for Meat and Carcass Traits Using Imputed Sequence Level Genotypes in Pooled F2-Designs in Pigs.

In order to gain insight into the genetic architecture of economically important traits in pigs and to derive suitable genetic markers to improve these traits in breeding programs, many studies have been conducted to map quantitative trait loci. Shortcomings of these studies were low mapping resolution, large confidence intervals for quantitative trait loci-positions and large linkage disequilibrium blocks. Here, we overcome these shortcomings by pooling four large F2 designs to produce smaller linkage disequilibrium blocks and by resequencing the founder generation at high coverage and the F1 generation at low coverage for subsequent imputation of the F2 generation to whole genome sequencing marker density. This lead to the discovery of more than 32 million variants, 8 million of which have not been previously reported. The pooling of the four F2 designs enabled us to perform a joint genome-wide association study, which lead to the identification of numerous significantly associated variant clusters on chromosomes 1, 2, 4, 7, 17 and 18 for the growth and carcass traits average daily gain, back fat thickness, meat fat ratio, and carcass length. We could not only confirm previously reported, but also discovered new quantitative trait loci. As a result, several new candidate genes are discussed, among them BMP2 (bone morphogenetic protein 2), which we recently discovered in a related study. Variant effect prediction revealed that 15 high impact variants for the traits back fat thickness, meat fat ratio and carcass length were among the statistically significantly associated variants.